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1 Introduction 



The propagation of quantum fields in curved spacetime is an interesting problem that 
is rich in surprises. For instance, it was shown by Drummond and Hathrell [1] (see also 
[2-7]) that vacuum polarization is sensitive to the curvature of spacetime and leads to 
a modification of the wave equation equation for photons in such a way that in the low 
frequency limit the phase velocity can actually be greater than c. 



The superluminal low frequency propagation of photons can be extracted from the 
effective action for QED in curved spacetime to linear order in the curvature: 
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The terms involving the curvature here are caused by vacuum polarization and the dots 
indicate that this effective action is the first term in a curvature/derivative expansion, 
so results deduced from it are valid only for low frequency propagation. At one-loop 
order, the Feynman diagrams that contribute are shown in figure 1 and it is graph (a) 
that is non-local and so sensitive to the background curavture. The effective action 
(1.1) has been renormalized and the one- loop divergence has been absorbed into the 
definition of the renormalized electric charge. 

The refractive index derived from (1.1) for a photon with wave- vector = cuk^, 
where u is the frequency, is 

n^j = Sij + — [ciSijR^^k^'k'' + CiR^^iujk^k'') , (1.2) 

where the indices i,j label the two spacelike polarisation directions. The constant 
coefficients ci,C2 are simply related to the known coefficients 02,03 in the effective 
action (see for example [10]). If we now introduce an affine coordinate u along the null 
geodesic with tangent vector k^, which we denote 7, we ca write the result above as 

nij{u) = 6ij + -^{ci6ijRuu{u) + C2Riuju{u)) . (1.3) 

This makes it clear that the Drummond-Hathrell result depends on the curvature com- 
ponents R\ju{u) and Ruuiu) evaluated along the null geodesic 7 describing the photon's 
classical trajectory. 

There are three scales in the problem, m and 9^, the latter being the scale of 
a typical element of the Riemann tensor of mass dimension [M"^]. The Drummond- 
Harthrell result based on the effective action in (1.1) is valid in the limit m? ^ $H and 
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Figure 1. The one-loop Feynman graphs that contribute to vacuum polarization. Only the 
first loop (a) is non-local in spacetime and so sensitive to the background curvature. 

cu^lH/m'^ ^ 1, while the description of the photon in terms of classical ray optics in 
curved space requires that cu^ ^ so that the wavelength is large compared with the 
length scale over which the curvature varies. 

The breakthrough made in a series of papers [8-11] was to extend the Drummond- 
Hathrell result by relaxing the requirement that u'^Dl/m'^ -C 1. The result is obtained 
by calculating the one-loop vacuum polarization due to the Feynman graphs illustrated 
in Figure 1 in the appropriate limits. The contribution from the first graph (a) is non- 
local in spacetime and therefore is sensitive to the background curvature. The effect 
of the non-local contribution corresponds to summing up an infinite set of terms in 
the effective action; however, this local expansion is generally only asymptotic and the 
complete result includes effects which are non-perturbative in the curvature. Note that 
the calculation is tractable because the two other limits, namely ^ 91 and ^ 91, 
are maintained, where the former means that we can still talk about the classical ray, 
or null geodesic 7, and the latter means that the result only depends on the metric 
in a tubular neighbourhood of 7. In fact this tubular neighbourhood can be made 
mathematically precise as the Penrose limit of the full metric associated to 7 [12] (see 
also [13, 14]). In terms of the so-called Brinkmann coordinates (m, f , z*) adapted to the 
geodesic 7, which is the curve v = = 0, the Penrose limit takes the form 

ds"^ = 2du dv - R\ju{u)z'z^ du^ + dz^ dz' . (1.4) 

The components R\ju{u) are both the components of Riemann tensor evaluated along 
7 and the non-vanishing components of the Riemann tensor of the Penrose limit itself. 
For later use, we also introduce the Rosen coordinates {u, v, x"), a = 1, 2, in which the 
Penrose limit metric (1.4) takes the form 

ds^ = 2du dv + Cab{u)dx'' dx^ . (1.5) 

The affine parameter u along 7 is common to both sets of coordinates. It turns out to 
be a great advantage to be able swop between the Brinkmann and Rosen coordinates 
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when necessary. This is achieved by introducing a zweibein for the transverse part of 
the Rosen metric 



Cabiu) = EUu)S'^E,h{u) , (1.6) 
specifically chosen so that the tensor 



% = ^E,^ (1.7) 



is symmetric. In the above Ei'^ is the inverse EiaEj"- = 5ij. In that case the relation 
between the two sets of coordinates is 

v = v + ^n,j{u)z'z^ , x" = E,\u)z' . (1.8) 

The expression for the 2x2 matrix of refractive indices in polarisation space is [10] 
n,,{u-u) = 5., + d^ai-OJ'.,[u; ^^^^^ _ , (1-9) 

with 

J-,, (m; z) = r 2e^^("'-") - A,, (m, u') v/A(m,m01 , (I-IO) 

In the above ^ is a Feynman parameter that is familiar in a one loop quantity. The cur- 
vature dependence is embedded in the "Van Vleck-Morette (VVM) matrix" Ajj («,«') 
and its determinant A(m, u') = det Aij{u, u'). This only depends on the curvature com- 
ponents R\ju{u) of the full metric evaluated along 7 and physically it encodes the tidal 
forces that are experienced along 7. In order to investigate these tidal forces, consider 
the geodesic equations for the transverse Brinkmann coordinates, or "Jacobi fields", 

^ + R\,^iu)z^ = 0. (1.11) 

The VVM matrix can then be defined in terms of the Jacobi fields as follows. If we 
write the solution for the Jacobi field z'^{u) in terms of some initial data at mV 

z\u) = ^,j{u,u')z^{u') + £^,j{u,u')z^{u') , (1.12) 



^Note that the transverse Brinkmann coordinates are raised and lowered with 6^ j and so one does 
not need to distinguish upper and lower indices. 
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then [10] 



Aij{u, u') = (m - u')£/jr\u, u') . (1.13) 

So the VVM matrix depends on the Jacobi fields that describe infinitesimal deforma- 
tions of geodesies that pass through u and u' on 7. The VVM matrix here, is of course, 
specifically defined in Brinkmann coordinates. However, the determinant A{u,u') is a 
bi-scalar quantity that can be written for arbitrary coordinates as 

A{u,u') = A(x(m),xK)) , A(x,x') = \ ^ det 'l^'^^'^j , (1.14) 

where a{x,x') is the geodesic interval between points x and x' (defined later in (3.8)). 
In Brinkmann coordinates it follows that 

A{u,u') = detAij{u,u') . (1.15) 

For later use, we quote an important identity [10] 

mj{u,u') + £^ik{u,u')nkj{u) = Eia{u)Ej''{u') . (1.16) 

Notice that A(m, u') is singular when u and u' are conjugate points on the geodesic 
7, that is when the geodesic through u and u' can be infinitesimally deformed. Hence, 
the integral in (1.10) is only defined with an appropriate prescription which involves 
deforming the contour into the upper half complex plane, as indicated. The integral in 
(1.10) is then convergent as long as the geometry becomes fiat in the asymptotic past 
u' —7- — oo. Notice that the result above for the refractive index is perfectly causal in 
that the refractive index at a point m G 7 only depends on the curvature along 7 in the 
past u' < u. 

The problem before us is to consider a similar philosophy for gravitons replacing 
photons, interacting with some massive scalar fields. We shall assume that the number 
of scalar fields is large enough so that the effect of matter loops is much larger 
than graviton loops, although in a more complete analysis graviton loops could also 
be considered. The analogue of the terms in the QED effective action (1.1) are those 
involving at most two powers of the curvature:^ 

S= ! d^x^\-^{R - 2A) + ciR^p^sR''^^' + c^Ra.pR''^ + csi?' + ■ ■ ■ 1 . (1.17) 
J LiovrG J 



^For example, see the monograph [16]. 
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This effective action is to be understood as being appropriately renormalized so that 
divergences are absorbed into the couphng G and the cosmological constant A. In 
d — A only two of the C(^H^) terms are independent as a consequence of the generalized 
Gauss-Bonnet theorem which implies that 



j d'x ^{R^p^sR'"'^' - ^Ra^R^' + R') 



[1.18) 



is a topological invariant whose addition to the action cannot affect the equations- 
of-motion. The three couplings q depend logarithmically on the mass m and the 
renormahzation scale /i. We will only need this dependence for the first two terms: 
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(1.19) 



The constants ctj are the arbitrary finite parts of the counter terms and manifest the 
fact that gravity coupled to matter fields is not a perturbatively renormalizable theory 
and so at each order in perturbation theory new couplings appear. 

In the next section, we investigate the effect of these terms of graviton propagation 
to find the gravitational analogue of the Drummond-Hathrell effect. 



2 The Gravitational Drummond-Hathrell Effect 

The classical propagation of gravitons in curved space is determined by expanding 
around a solution to Einstein's equations to linear order g^,^ g^,^ + h^^. As usual, 
one needs to fix a gauge in order to uncover the physical degree-of-freedom and to this 
end we impose the conventional transverse traceless gauge: 

V^h^"" = , g^^h^^ = . (2.1) 



In a general background spacetime, it is not possible to solve for the graviton modes 
exactly. However, just as in the QED case, we will work in the WKB approximation, 
which is valid when the frequency is much greater than the scale over which the cur- 
vature varies, ^ •[H. In this case, we can write the metric perturbation in the 
form 

Ku{x) = £^.(a;)e^®(-) , (2.2) 
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where the eikonal phase 9 is and the polarization tensor e^^ is 0{uj^). Substi- 

tuting into the equation-of-motion, and expanding in powers of l/w, the leading order 
term yields the eikonal equation for the phase: 

ae ■ ae = ; (2.3) 

which implies that the gradient = d^Q is a null vector field. This vector field defines 
a null congruence, that is a family of null geodesies whose tangent vectors are identified 
with the vector field k'^. This vector can also be identified with the 4-momentum of 
photons: in this sense the eikonal approximation is the limit of classical ray optics. 

We can then pick out a particular null geodesic, or classcal ray, in the congruence. 
Our goal is to find how this particular ray 7 is affected by the 0{D\'^) terms in the 
action. There are a set of coordinates {u, v, x"") which are specifically adapted to the 
null congruence: these will become the Rosen coordinates of the Penrose limit (1.5). 
In these adapted coordinates the geodesies are simply lines of constant {v, x"') and u is 
the affine parameter. We will choose the ray 7 to be the one with x"' = v = 0. In these 
coordinates the eikonal phase is simply, 

e = ujv . (2.4) 

As explained in ref. [13], the full metric Q/j^j, can always be brought into the form 

ds'^ = 2du dv + C(u, {), x"')dv^ + 2C„(u, {), x^)dx'' dv + Cab{u, 0, x'')dx'' dx^ . (2.5) 

In order to find the effect of the 0{y^) terms on the gravitation propagation we 
can then expand in the transverse coordinates (5, x"") around 7. In fact, the terms that 
we need in order to fully describe the effect is precisely the Penrose limit of the metric 
around the null geodesic 7. This correspond to expanding in {v,x°') and keeping the 
terms in (1.5). The adapted coordinates are then identified with the Rosen coordinates 
of the plane-wave metric. The Penrose limit has precisely the information needed to 
describe the behaviour of geodesies in the neighbourhood of 7 to linear order and we 
will see precisely encodes the information needed to investigate the effect of curvature 
on graviton propagation in the eikonal limit. 

Once we have approximated the metric with Penrose limit, we can solve for the 
graviton modes exactly, these modes then correspond to the graviton modes of the full 
metric to leading order in the WKB approximation. There is a subtlety here in that 
to capture the leading order approximation involves keeping some of the components 
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of the polarization tensor e^i, which are order but also other components which are 
order u^^. This becomes clear when one solves for the exact complexified graviton 
modes in the plane-wave background: 

_ r'E.'(u)E/(u) ^ _ . P«2s(2Le- . (2.6) 

with the remaining components vanishing. The two independent polarization tensors 
Ps, s = 1,2, can be chosen as 

Notice that, as alluded to above, at leading order in the WKB approximation, we need 
to keep the component h^^ even though it goes like uj~^. 

In order to describe the quantum corrections to graviton propagation, one takes 
the graviton modes off shell by modifying the eikonal phase via a u dependent factor: 

&-^oo{v + e,Au)) . (2.8) 

The corrected phase is a 2 x 2 matrix in the two-dimensional polarization space s = 1, 2. 
In general, 6ss'{u) will have both a real and an imaginary part. The real part describes 
a local correction to the phase velocity of the graviton along 7 and the imaginary part 
describes a local correction to the amplitude of the graviton along 7. Since Oss'iu) is 
perturbatively small, we can think of the result in terms of a matrix of refractive indices 

n..K«;^) = 5..' + 2%^^. (2.9) 

du 



The one-loop quantum corrections to the equations-of-motion that are linear in the 
curvature come from the last three terms in (1.17). It turns out that only the first two 
terms of order 0{d^) in (1.17) contribute to graviton propagation. The Gauss-Bonnet 
argument can then be used to reduce this to a single term involving the square of the 
Ricci tensor. So for a plane wave metric the renormalized effective action to this order 
takes the form 



S 
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(2.10) 



for a finite quantity A which is an independent dimensionless coupling constant. Sub- 
stituting the mode (2.6) with modified eikonal phase (2.8) into the equation-of-motion 
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that follows from the above action at leading order in the eikonal approximation yields 
a simple equation of the form 



where we have defined the effective Newton constant Gefr = GN.'^ Just as in the QED 
case, the result can be stated as a curvature induced contribution to the refractive 
index for graviton propagation along the geodesic 7, 

risAy^) = (1 + ^^-(")) ' (2-12) 

to 0{Gcs)- The result is actually simpler than the QED case in (1.3) since the result 
is diagonal in the 2-dimensional polarization space. So the effective of curvature, in 
the limit ^ SH, is to induce a frequency independent shift in the refractive index. 
Notice that this effect depends on the component of the Ricci tensor Ruu{u) of the 
Penrose limit which is equal to the same components of the Ricci tensor of the full 
metric evaluated along 7. 



3 One-Loop Vacuum Polarization 

The aim of this section is to compute the effect of one-loop vacuum polarization of a 
massive scalar field on graviton propagation in the eikonal limit cu^ ^ £H and in the 
limit that ^ fK. This latter limit, means that the quantum fluctuations around 
the classical geodesic 7 are small and so allows us to approximate the metric with 
the Penrose limit around the geodesic 7. This leaves a non-trivial dependence on the 
dimensionless ratio cj^D^/m'*. In particular, when we talk about low and high frequency 
we always mean the limits: 

low: — -T- < 1 > high: — - > 1 . (3.1) 

but subject to u"^ ^ 

At the one-loop level there are two graphs that contribute to the quantum equation- 
of-motion of the graviton field. The first is the important one from our point-of-view 
because it is non-local and therefore sensitive to the curvature. The second is local 
and we will be able to account for it by subtracting the contribution of the first in the 

■^We are assuming that » 1 so that graviton loops are suppressed relative to matter loops. 
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flat space limit. So turning to the first graph. It involves the 2-point function of the 
energy-momentum tensor coupled to the external gravitons. Since we are interested in 
the quantum corrected equation-of-motion, the relevant 2-point function is the retarded 
one 



Ilf,uap{x,x') = local + (0|T^^(a;)T^p(x')|0)i.et • (3.2) 

The fact that the retarded 2-point function is needed ensures that the quantum equation- 
of-motion is causal, since it vanishes for x' lying outside the backward lightcone of x. 
In the above the energy-momentum tensor of the scalar field is 

T,, = d,<l)d,<P - \g^u{dP<pd,(j) + m'<f)') . (3.3) 

The linearized quantum-corrected equation-of-motion for the graviton modes takes the 
form 

□V(x) = 87rGeff j d^x' ^-g{x') n^,,,(x, x')h''P{x') , (3.4) 

with, on the left-hand side, the appropriate Laplacian for gravitons in a non-trivial 
background. The strategy is to solve this order-by-order in the coupling Ges- At 
leading order, one takes the classical graviton mode (2.6) with the corrected eikonal 
phase (2.8) with 6 of order Geg one the left-hand side and on the right-hand side the 
classical graviton mode. This leads to the equation for the one-loop 0{Gqs) correction 
to the phase: 

^(^ifM = local + f c/V ^^^)hl^'''{x)d^d'Mx,x') d,d'fi{x,x')h1f{x') . 

(3.5) 



The insertions of the classical graviton modes involve, from (2.6) 



A/detE^ 



(3.6) 



Note only the first term in the energy-momentum tensors contribute due to the gauge 
condition (2.1). In the above, G{x,x') are the scalar propagators in the plane-wave 
background. These can be taken as Feynman propagators even though this does not 
give the retarded contribution to the 2-point function of the stress tensor required in 
(3.2). The reason why we can make this simplification is that we actually want the 
convolution of the 2-point function with the classical graviton mode as in (3.4). The 
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latter is a positive frequency solution with respect to the null coordinate u and this has 
the effect of picking out the retarded part of the 2-point function of the stress tensor.^ 

The scalar Feynman propagator in a general background spacetime can be written 
in the heat-kernel or "proper-time" formalism as 

,e-^™ ^+2T'^(-'-')(](a;, x'\T) , (3.7) 

subject to the usual rn? — ie prescription. Here, a{x, x') is the geodesic interval 

between the points x and x': 

1 /•! 

a{x,x') = - dr g^y{x)x^x'' , (3.8) 
^ Jo 

where x^ = x^(r) is the geodesic joining x = x{0) and x' = x{l). The factor A{x,x') 
is the famous Van Vleck-Morette determinant defined previously in (1.14). 

The expression (3.7) has a nice interpretation in the worldline formalism, in which 
the propagator between two points x and x' is determined by a sum over worldlines 
x'^{t) that connect x = x{0) and x' = x{T) weighted by expiSfs] with the action 

1 

S[x] = -m^T + - dr g^,{x)x''x'' . (3.9) 
4 Jo 

Here, T is the worldline length of the loop which is an auxiliary parameter that must 
be integrated over. The expression (3.7) corresponds to the expansion of the resulting 
functional integral around the stationary phase solution, which is simply the classical 
geodesic that joins x and x'. In particular, the classical geodesic has an action S[x] = 
a{x,x') /2T — m?T giving the exponential terms in (3.7). The VVM determinant comes 
from integrating over the fluctuations around the geodesic to Gaussian order while the 
term VL{x^x'\T) = 1 + Yl^= 

')T^ encodes all the higher non-linear corrections. 
These terms are effectively an expansion in ^H/m^, so the form for the propagator 
is useful in the limit of weak curvature compared with the Compton wavelength of 
the scalar particle. Of course, this is precisely the limit we are working in here and 
this goes hand-in-hand with the fact that in the plane wave limit this factor is trivial 
Vt{x,x'\T) = 1. 

The weak curvature limit 9^ ^ leads to a considerable simplification as we now 
explain. Denoting the worldline lengths of each of the scalar propagators as Ti and T2, 

"^This important point is explained in detail in [10]. 
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the terms in the exponent of the integrals in (3.5) are 



exp 



^rn^T + 1 (^1 + 1) x') + lujv'] • (3.10) 



For later use, we find it convenient to change variables from Ti and T2 to T = Ti + T2 
and ^ = Ti/T, so < ^ < 1. The parameter ^ is identified as a conventional Feynman 
parameter. Expressed the other way 

Ti = Te, T2 = T{l-i). (3.11) 

The Jacobian is 

ciTi dT2 r dT ^2) 







Ti Jo T^Jo m-0? ■ 



In the limit 9^ ^ the integral over x' is dominated by a stationary phase 
determined by extremizing the exponent (3.10) with respect to x': 

^^^^i-^a;a(a:, x') + ud'^v' = . (3.13) 

Since d''^a{x,x') is the tangent vector at x' of the geodesic passing through x' and x, 
the stationary phase solution corresponds to a geodesic with tangent vector a d'^^v'. 
This means that x and x' must lie on a null geodesic. If we choose x to be the point 
(m, 0, 0, 0) then x' must have Rosen coordinates («', 0, 0, 0) and so x' and x lie on the 
geodesic 7 and 

di,'a{x,x') = u — u ; (3-14) 

so the v' component of (3.13) becomes 

u' ~ u , ^ 

+ w = (3.15) 



2re(i-0 

and hence 



u =u-2ujTi{l-Cj ■ (3.16) 



In the equivalent worldline picture, the stationary phase solution which dominates 
in the limit ^ m? describes a situation where the incoming photon propgating along 
7 decays to an electron positron pair at the point u' = u — 2ujT^{1 — ^) which then 
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propagate along the null geodesic 7 to the point u and then combine into the photon 
again. This was a key step in the derivation of the refractive index in the worldline 
formalism in [8, 9]. 

The plane-wave limit of the full metric provides precisely the data that is needed 
to perform the saddle-point approximation to leading order. In other words, for a 
plane wave background the problem is exact at Gaussian order. For the remainder of 
the calculation, we now simply assume a plane-wave metric knowing that the result 
will apply to the full metric when £R; ^ m^. The meat of the calculation involves 
performing the Gaussian integral over x' where the exponential factor has the form 
(3.10). In the plane-wave background the geodesic interval has a simple form quadratic 
in the transverse coordinates 

a(x,x') = {u~u'){v-v') + \/\ab{u,u'){x-x'Y{x-x'f , (3.17) 

where 



A„,,(u,m') = (u - u') 



n -1 



C-\u")du" 



(3.18) 



ab 



is the transverse part of the VVM matrix in Rosen coordinates. It is related to the 
same quantity in the Brinkmann coordinates by the zweibein 

^ij{u,u') = Ei\u)Ej\u')^ab{u,u') . (3.19) 

The integrals over v' and u' are trivial, the former leading to a delta function constraint 

which then saturates the integral over the latter. The condition on u' is precisely the 
saddle-point condition (3.16). 

It remains to perform the Gaussian integrals over the two transverse coordinates 
x'"". The Gaussian integrals that we need are of the form^ 



r(l) 



J A/det Aafe(n,u') 



^Where the integrals are rendered convergent by the prescription T ^ T — zO^ and we define 
x' ■ A(m, u') ■ x' = x"'Aabiu, u')x"'. 
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for (3.5) with indices {/i, u, a, p} = {v, v, v, v}, along with 

(3.22) 

for (3.5) with indices {fi,h',a,p} = {i,j,v,v} or {v,v,k,l}, and 

r(3) ^ / J2 f d'^ ^x'-A{u,u'yx' ^^f\_^^ x'-A{u,u'yx' 
abed J dx'^dx'^ dx'^dx"^ 

= (A..(^, u')AM, W) (3-23) 

T^det /\ah{u,u') 

for (3.5) with indices {/x, z/, a, p} = {i,j,k,l}. 

The contribution for indices {i,j,k,l} involves the quantity 

E,%u)E,\u)lifju, u')E,\u')Ei\u') , (3.24) 
which, using (3.19) is proportional to 

Ej''{u)E,,a{u')A,j,{u,u')Ayi{u,u')Ek>biu)Ek\u') 

(3.25) 

+ Aik{u, u')Aji{u, u') + Aii{u, u')Ajk{u, u') . 

The contribution form the other terms involving indices {;,{)}, {v^v^k^l} and 
{{), -0, £>, £>} is proportional to 

ni,{u)l'S{u,u')Ek'{u')Ei\u') + E,\u)E^\u)I^;S{uy)^ki{u') 

i(u-u'V .^ (3-26) 

which, using (3.16), (3.21) and (3.22), is proportional to 

nij{u)£/i^,\u',u)Ek4u)Ek\u') + Ej\u)Ej,a{u')£^yl{u\u)nki{u') + nij{u)nki{u') . 

(3.27) 

Using the identity (1.16), and being careful with the relative normalization, these terms 
combine nicely with the first term in (3.25), to give a net contribution involving the 
tensor^ 

1 ... .^T. .^ X (3-28) 



{u — u') 



'Here, we use the symmetry properties si^ij{u,u') — —£/ji{u\u) and SSij{u,u') — Mji{u',u). 
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In particular, if we had not included the h^"" term in (2.6), we would not have been 
able to write the final result in terms of the matrices ^ and ^ and, for example, 
the result for a symmetric plane wave to be discussed below would not have involved 
a function of the difference u — u' which on physical grounds we expect since this 
particular plane wave admits a Killing vector du- But with the contribution from K"^ , 
the final expression depends only on the matrices £^ij{u, u') and ^^ij{u, u') or Aij{u, u') 
that are associated to the Jacobi fields. The final tensorial nature of the result can be 
summarized by the quantity 

dcf. 1 



AgL(^,^') 



16 



(A(u, u')^^(m, u'))ij{BS^ {u, u')A(m, u'))ki 



+ Aife(M, u')Aji{u, u') + Aii{u, u')Ajk{u, u') 
Collecting all the remaining factors together gives the result 

local 



(3.29) 



du 



-im2rA{2) 
ss 



u, u 



(3.30) 



u' =u-2uiT^{l-e,) 



where 



A 



S(n,nO=P^Ag),(«,«')P^^ 



(3.3i; 



The local contribution comes from graph (b) in Figure 1. We can account for this 
contribution by subtracting the fiat space limit of the non-local expression. This gives 

foo— iO+ 



du 



eff 



167ra;^ 

X 



d^e-" 



dT 

AfJ,{u,u')^/A(^' 



(3.32) 



J u'=u-2ujTi{l-i) 

The contour for the T integral is taken just below to real axis in order to avoid the 
singularities of the VVM determinant when u and u' are conjugate points. 

The resulting integral is still divergent as we can see by expanding the integrand 
in powers of uj. This is achieved by expanding in powers of m — u': 



Afj {u, u') ^jA{u,u') = (4 + Ruu{u) {u - u'f) + 0{u- u'f 



This means that divergent contribution is of the form 



du 



off 



div. 



1207r 



Ruu{u)6s 



dT 



-im^T 



(3.33) 



(3.34) 
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We can regularize by first rotating to Euclidean proper time T — —iT and then 
introducing a cut off 5 



dOss'iu) 



du 



G. 



cff 



div. 1207r 



7i5 + log(m^5) Ruu{u)6, 



(3.35) 



We can subtract the divergence and introduce an arbitrary finite piece in the form of 
the coupling A. The regularized expression is then 



dT 



d9 



i-m?T 



X 



Afj(u,M') yAl^M^ - (4 + Ruu{u){u - u'f)5ss' 



(3.36) 



J «'=u-2wTg(l-^) 

where the integrand is now regular &t u' = u. In the low frequency limit, we have 

d9.,.,'(u) XGcs „ ... , m'^Gcs f u;V^\ . . 



du 



which matches the form extracted from the effective action (2.11) precisely. 

We can express the result above as a matrix of refractive indices and change vari- 
ables from T to u' to write the result in the same form as the QED result (1.9): 



nss'iu;ij) = 6s 



G 



cff 



27r 



^ d^m-OfJ'ss'lu; 



m 



(3.38) 



where 



Tss'iu; z) = -\Ruu{u)5, 
du' 



+ 



oo+i' 



0+ {u - u'f 



iz{u' —u) 



Afj(n, u')^/\{uy)- (4 + Ruu{u){u - u' f))5ss' 



(3.39) 



This result is the gravitational analogue of the formula for the refractive index of the 
photon in (1.10). 



4 Properties of the Gravitational Refractive Index 



In this section, we investigate some of the phenomenology of the gravitational refractive 
index. In order to have a concrete example to hand, we first look at the simplest kind 
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of plane-wave spaces for which the Riemann tensor components R^uju are constant 
along the geodesic. These are the symmetric plane waves, or Cahen-Wallach spaces 
[17]. To simplify things further, we take the case where the two eigenvalues of the 
Riemann tensor are equal, so up to a choice of coordinates we can write R\ju = (J^Sij, 
for constant a. In this case, the VVM matrix, and related tensor A*^^\ are 



A,, = —-. '—5,j , A J = '-5,s' . 4.1 

sm.a[u- u') a(u-u') 

In this case, everything is diagonal in polarization indices and so we will not show them 
where possible. Changing variable to t = u — u\ we have, from (3.39), 



oo— iO+ 
2 I / „-izt 



J^iz) = -2Xa' + dte 



-sin^crt t J 2) 

-2Xa^ -Aza + ^ + 2{z^ - a^) [log j + ^(^^— j 



where il){x) is the di-gamma function. Notice that the integrand is a function of the 
difference t = u ~ u' which is a special feature of a symmetric plane wave background 
which is translationally invariant under u ^ u + a. 

The first issue is to investigate is the behaviour of the refractive index at high and 
low frequencies, or more precisely large and small values of the dimensionless ration 
bJG jwt} . The low frequency behaviour is obtained by expanding (4.2) in powers of 2;"^: 

"(^) = ' + ^1^ + 472^ + 945945^ ' ^'"'^ 

So the low frequency limit of the refractive index is independent of w and depends 
on the coupling constant A. Consequently, the low frequency phase velocity is not 
necessarily equal to c. As has been discussed in detail in [7], it is possible that n(0) < 1 
since this does not imply a violation of causality because low frequency waves cannot 
be used to send information. 

Causality actually depends on the high frequency limit of the refractive index, since 
this is what governs the propagation of a sharp wavefront. For our plane- wave example, 
this can be extracted by taking the z — )■ limit of (4.2), 

T{z) 2a^ log z + 0{z^) . (4.4) 

Therefore, the refractive index behaves as 



n[uj) = l + — — loga; + --- , (4.5) 
oOtt 
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for large u. This indicates that at some high frequency the perturbative expansion is 
breaking down and — at the very least — one would need to sum up an infinite class of 
higher terms in perturbation theory in order to get a sensible result. The result above, 
for the high frequency behaviour is in fact universal, since the origin of the logarithm 
is in the divergence at the upper limit of integration of the third term in the integrand 
(3.39) as z —¥ 0. The general result for the high frequency limit, is therefore simply 

nss'iu;u) '^=°° 6 ss' + -^7^ Ruuiu) log u 6 ss' -\ • (4-6) 

Interestingly, the logarithmic behaviour is increasing as long as the null energy condition 
is satisfied since this requires that Ruu > 0. 

For our simple example of a symmetric plane- wave space, the refractive index is 
real, however, more generally the refractive index has both a real and an imaginary 
part. The real part represents a change in the phase velocity along the ray 7 while the 
imaginary part represents the dissipative effect of propagation caused by decay into 
pairs of the scalar particle, although there is nothing to prevent an amplification, as we 
explain below. For instance, if we expand the result (3.38) in powers of the curvature 
by taking the next term in the series (3.33) 

= (4 + Ruuiu)iu - n')' + ld^Ruuiu)iu - v!f)b,, + 0{u - , 

then one finds to order fH^/^,^ that there is a contribution to the imaginary part of the 
refractive index involving the derivative of the Ricci tensor: 

Im nss'{u\ u) = 2^^^°Q^2 ' duRuu{u)5ss' H • (4-8) 

This represents a local quantum contribution to the amplitude of the graviton involving 
the Ricci tensor of the form 



exp 



27r ■ 140m2 



Ru,,iu) 



(4.9) 



as the graviton propagates. This change in the amplitude is caused by the effect of 
the local curvature on the wavefunction renormalization of the graviton [18]. In order 
to understand this effect we have to remind ourselves that a physical graviton is a 
renormalized composite object that consists of a bare graviton plus a cloud of virtual 



'^Note that a derivative of a component of the Riemann or Ricci tensors is order ^{■^/'^ in our 
convention. 
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scalar pairs. In flat space, there is a delicate balance between the bare graviton and the 
virtual cloud of scalars, but in curved space this balance is perturbed by the tidal forces 
acting on the composite particle induced by the curvature of spacetime. The quantity 
r = 2a; Im n{u]u) represents the instantaneous rate of production of scalar particle 
pairs from the gravitons. It is important to realize that this can either be positive or 
negative as a function of u, since the balance between bare particle and virtual cloud 
can be upset either way. In particular, in certain backgrounds the graviton can actually 
amplify as the virtual cloud becomes depleted. This does not break unitarity or violate 
the optical theorem because the graviton is already a renormalized particle. One way 
to investigate this, is to imagine turning on the coupling at finite affine time u = uq. 
In this way, one would see the graviton undergo wavefunction renormalization even in 
fiat space. However, in a gravitational theory, as in QED, this wavefunction renro- 
malization is a UV sensitive process that actually diverges as the cut-off is removed.® 
So the dressing which creates the renormalized graviton is an infinite process and the 
implication of this is that the graviton may become infinitely "undressed" which means 
that the amplitude can increase without bound as the graviton becomes undressed by 
curvature induced tidal forces. More generally the imaginary part can be written in 
nice way by changing variables io t = u — u' and extending the contour to run from 
— oo to +00: 

Im Uss'iu; io) 

Jo J-oo-iO+ ^ 

A useful consistency check of the result (3.39) can be obtained by taking the ex- 
ternal graviton off-shell which is possible in a plane-wave background since it can be 
done simply by modifying the eikonal phase in (2.4) as follows 

Q = ujv — >ujv^ , (4.11) 

2uj 

where M is the off-shell mass of the graviton. This extra phase can be carried through 
the one-loop computation of the graviton propagation. This extra phase only affects the 
non-local contribution from Feynman graph (a) in Figure 1 and leads to a modification 
of the integrand in the refractive index formula in (3.39) of the form^ 

du' 



L 



Afj(M,M')v/AKwO • (4.12) 



^This effect is deseribed in detail in [15]. 

^Note that the subtractions in (3.39) which become the local contribution to the vacuum polariza- 



tion are not modified and they are independent of M. 
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The expression for the imaginary part of the refractive index in (4.10) is similarly 
modified. Now when M > 2m, which corresponds to the off-shell graviton being above 
threshold for decay into pairs of the scalar particle, the imaginary part of the refractive 
index can be evaluated by completing the contour in (4.10) in the upper-half t-plane.^° 
In particular, there is a contribution from the pole at t = 0, which is independent of 
the curvature, of 

Inin(u;) = !I^^ C d^e{M^i{l - ^) - m^) = "^^^^^1 - Am'' /M^ , (4.13) 

diagonal in polarization indices, which simply represents the exponential decrease in 
the off-shell graviton amplitude in flat space due to decay into pairs of the scalars above 
threshold. The decay rate is given by 2uj Im n in agreement with a conventional fiat 
space calculation of the rate. 

There is another way that an imaginary part of the refractive index can be produced 
which is non-perturbative in the curvature. For a symmetric plane wave where at least 
one of the eigenvalues of the constant BJ'uju is negative there is a positive contribution 
to Imn(ci;), so a decaying amplitude, which is non-perturbative in the curvature. The 
simplest case to consider is BJ'uju = — c^^ij- Such a background violates the null energy 
condition, but provides a useful example where the imaginary part of the refractive 
index can be calculated exactly. In this case, \m.n{uj) = to all orders in the curvature 
expansion as can be seen from (4.3) by taking a — ?■ ia. In this case, the VVM matrix 
and related tensor A*^^) are obtained by substituting cr — )■ io" in (4.1) 

smh a{u - u') smh^ o-{u - u') 

The integral in (4.10) can be computed by deforming the contour and picking up the 
contributions from the poles on the negative imaginary axis at t = —irni /a, n = 1, 2, . . .. 
This gives directly 

which is, indeed, non-perturbative in the curvature. In this case the result is indepen- 
dent of u and so it implies a constant rate of production of pairs of the scalar particles. 
Such a process is kinematically disallowed in fiat space, but there is no such kinematic 
constraint in curved space. Although this result has been calculated in a background 



-^°More precisely only for the values of ^ for which M^^(l — ^) > m^. 
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that violates the null energy condition it is generic to plane waves with at least one 
negative eigenvalue. 

To conclude it is interesting to compare the "phenomenology" of the gravitational 
refractive index with the photon refractive index described in detail in [10, 11]. The 
detailed shape of the real and imaginary parts of the refractive index depends sensitively 
on the curvature in both cases. The obvious difference is that the photon refractive 
index approaches unity at high frequencies like l/w, to compare with the logarithmically 
diverging behaviour of the graviton case. This is reflection of the fact that QED is 
better behaved in the UV compared with quantum corrections to general relativity 
from matter fluctuations. In particular, the one-loop QED calculation is entirely self- 
consistent whereas the one-loop calculation in the graviton case predicts its own demise 
at sufficient high frequency. 
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